The misleading nature of the leaky box models in cosmic ray physics 
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Abstract: Many experimental results around and above the energies where the solar modulation affects 
cosmic ion fluxes were quantified, conceptualized and debated using leaky box models. These models 
exploit the notion of equilibrium between creation and destruction processes of cosmic ions in an undif- 
ferentiated arbitrary volume representing the Galaxy, ignoring the galactic magnetic field, the size of the 
Galaxy, the position of the solar cavity, the spatial distribution of the sources, the space variation of the 
interstellar matter and other pertinent observations. Progress in the measurements of the quoted observa- 
tional parameters rendes obsolete the use of the leaky box models. Specific examples substantiating the 
inadequacy of the leaky box models are analyzed such as the conversion of the boron-to-carbon flux ratio 
into grammage and the residence times of cosmic ions in the Galaxy. The unphysical and misleading 
nature of the leaky box models is ascertained and illustrated at very high energy. 
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Introduction 

Many experimental results on the primary cos- 
mic radiation at very low energy gathered in the 
last five decades by balloon-borne instruments and 
satellite experiments have been expressed using 
simple theories of galactic cosmic rays referred to 
as Leaky Box Models. Two basic physical quanti- 
ties, namely the gas column encountered by Galac- 
tic cosmic rays (grammage) and the related resi- 
dence time (T), were believed to correctly inter- 
pret a number of measurements. Ritual fittings of 
the grammage above an adjustable energy E$ (or 
the rigidity Rq) based on numerous measurements 
have the form: 



.9 = Gr 



-6 



(1) 



where G is a constant grammage at the rigidity Rq, 
8 a constant and r = R/Rq, with R the ion rigid- 
ity. Numerical values in classical interpolations 
are, for example: G=10.8 g/cm 2 , <5=0.6 and i?o=4 
GV [19], or G=24.0 g/cm 2 , <5=0.65 and i? =5.5 
GV [12]. 

Grammage is converted into residence time by: T 
= glpv where v is the ion velocity and p the mat- 
ter density. The extrapolation at very high energy 
of the residence time caused fictitious problems, 
as pointed out by Hillas [4, 20]. Figure 1 delib- 
erately shows on a linear scale in energy the B/C 
flux ratio in the energy band 10 - 200 GeV/u indi- 
cating no compelling empirical evidence for a de- 
crease of the grammage of the form r~ s which, on 



the contrary, is well established in the region 1 — 4 
GeV/u. The recent measurements of the B/C flux 
ratio at 700 GeV/u by the Runjob Collaboration 
[18], along with the data shown in figure 1, exclude 
the functional form [19] above 20 GeV/u. The data 
in figure 1 and the B/C flux ratio [5] vividly tes- 
tify to the inadequacy of the grammage fitting via 
r~ 6 to physical reality in an area believed in past 
decades to be the realm of the Leaky Box Models. 



The unjustified form of the cosmic-ray 
source power versus energy 

In a suitable Galactic container, at a given energy 
E, two processes are at work to maintain cosmic- 
ray intensity at a constant value: escape from the 
container and nuclear interactions inside. The sim- 
plified equation describing these processes are: 

dQ/dE = k{dN/dE)(l/X + l/f) (2) 

where dQ/dE 1 is the source power, dN/dE is the 
differential intensity in a given location of the disc, 
A is the nuclear collision length, / is the average 
escape length and k a suitable normalization con- 
stant. It is an established result that the differential 
intensity at Earth, dN/dE, measured in many en- 
ergy bands, obeys a power law i.e. dN/dE = aE~^ 
where 7 is the spectral index and a a constant. 

Whenever one of the two terms (1/A) and (1//) 
dominates, the equation (1) simplifies further. At 
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Figure 1: Measurements of the boron-to-carbon 
flux ratios above 10 GeV/u up to the energy of 200 
GeV/u displayed in a linear scale of energy. No de- 
creasing trend of the B/C flux ratio is evident. Data 
below 10 GeV/u may be found elsewhere [22]. 

sufficient high energy, especially for light ions, the 
term 1/A is negligible compared to 1//. Using the 
ad hoc hypothesis that the escape time has the form 
E~ 5 it follows that dQ/dE is proportional to E~ s 
where s= 7 - S. Therefore from the assumption 
E~ 5 (disproved by the B/C data above 10 GeV/u) 
the unknown form of the source power versus en- 
ergy of the cosmic radiation automatically trans- 
forms into a power law with a constant index s. 
The splitting of 7 in two constant parts, s and 5, in 
a large energy band remains unjustified taking into 
account, not only the B/C flux ratio at 160 and 700 
GeV/u measured by the Runjob experiment [18], 
but also the classical data shown in figure 1 . 

The unphysical grammage extrapolated 
at high energy 

From simple elaborations of the equation (1) it 
follows that the grammage, g, is related to the 
secondary-to-primary flux ratio of cosmic rays (for 
example, B/C but also 3 He/ 4 He, subFe/Fe, Nitro- 



Figure 2: Interpolated grammage via G r~ s (thin 
lines a and b) for two sets of constants mentioned 
in the Introduction. It is evident that the slope S 
tuned to the B/C flux ratio in the range 1-5 GeV/u 
is not in accord with the experimental data above 
10 GeV/u. The computed grammage (c) is a tiny 
fraction of that displayed in figure 3. 

gen/Oxygen and others) by the equation: 



where the N s /N p is the flux ratio observable by 
experiments at a given energy, m the mean mass 
of the interstellar atom, <r p is the inelastic nu- 
clear cross sections of the primary with the inter- 
stellar matter and a sp secondary production rate. 
Equation (3) applies to a single parent nucleus 
generating a unique secondary while primary-to- 
secondary ratios resulting from many tributaries 
have more articulated formulae. 

In the limited energy range, 1 - 5 GeV/u, there is 
empirical evidence that N S IN P ratio (i.e. the B/C 
ratio) is decreasing with energy. Consequently, 
since <r sp and a p are smooth functions, the gram- 
mage has a decreasing trend with energy, often 
interpolated by Gr~ 5 with a constant S. On the 
contrary, in the energy band 10 - 200 GeV/u, 
the classical measurements of the B/C flux ratio 
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Figure 3: The average distance of the sources from 
the Earth in the energy range 10 10 - 10 18 eV and 
the related grammage of galactic cosmic rays (He 
and Fe) intercepting the Earth, expressed in g/cm 2 
in the same vertical axis. This calculation, unlike 
those based on Leaky Box Models, exploits numer- 
ous astronomical and radioastronomical observa- 
tions (see Section 2 in [2]). 

[16, 11, 8, 15, 9, 7, 14] do not conform to the same 
S extracted at lower energies, as apparent from the 
data shown in figure 2. The functional form r~ s 
fails just above 20 GeV/u, exhibiting a large dis- 
crepancy at 700 GeV/u. The B/C flux ratios mea- 
sured by the Runjob experiment [18] reconfirm the 
failure of the formula (1) with a unique 5 in the 
interval 1 - 700 GeV/u. 

It is surprising that most of the data in the region 10 
- 100 GeV/u shown in figure 1 are simply omitted 
in the comparison with the formula g = Gr~ 5 (see, 
for example, fig.l [13] or fig.l [10]) or with direct 
calculations (fig.2 [5]). 

The grammage, the gas column encountered by 
Galactic cosmic rays, g, can be calculated using 
directly the trajectory length Lrj, via g = m n 
Lr>, where n the number density of the interstellar 
atom. Figure 3 gives the He and Fe grammage ver- 
sus energy in the interval 10 10 — 10 18 eV. Details of 
the calculation are given elsewhere (13). The com- 



puted He grammage diminishes with a gentle, in- 
creasing slope up to 10 15 eV beyond which a steep 
descend sets on, finally attaining a minimum value, 
at 10 18 eV (rectilinear propagation regime). The 
dependence of the He and Fe grammage versus 
energy is explained elsewhere [1, 2] by the notion 
of a Galactic basin. 

The grammage computed by the direct formula, g 
= m tih Lrj, is certainly different, by definition, 
from that evaluated by the transport equation used 
in Leaky Box Models, which relates the grammage 
to the B/C flux ratio. For instance, the spatial dis- 
tribution of the secondaries should differ from that 
of the primaries and this difference is not incorpo- 
rated in Leaky Box Models. The carbon grammage 
versus energy evaluated by the trajectory length, 
via g = m nn in, is shown in figure 2 (thick line) 
arbitrarily normalized at 50 GeV/u to the B/C ratio 
of 0.087. The slope of the grammage versus energy 
in the band 50 - 700 GeV/u is consistent with the 
data. Note that at lower energy, below 10 GeV/u, 
the grammage for stable 7 Be has been evaluated 
with the same method resulting in a much steeper 
slope [21]. 

Additional calculations [3] of the grammage for 
the stable 9 Be result in a grammage decrease of 30 
per cent in the range 1-10 GeV/u both for the spi- 
ral and the circular galactic magnetic fields. Note 
that the mass and the charge of 9 Be are similar to 
those of B and C, and consequently, the grammage 
profiles with energy are analogous. The measured 
B/C flux ratio decreases from 0.32 at 1 GeV/u to 
0.20 at 10 GeV/u, exactly the same decrease of the 
computed 9 Be grammage (see fig.2 [3]). 

The unphysical walls reflecting back 
cosmic rays in the disc 

The ad hoc assumption that the residence time of 
cosmic rays is compatible with a single value (for 
example, 15 x 10 6 years) is unphysical because cos- 
mic rays originated in the Bulge resides longer that 
those populating the disc periphery, at 15 kpc from 
the galactic center. The volume where cosmic rays 
propagate is not specified and therefore the resi- 
dence volume is undefined. The physical motion 
(migration, diffusion, convection, trapping or com- 
binations of these classes of displacements) of cos- 
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mic rays is not defined. Since the magnetic field 
does not exist in Leaky Box Models and the ion 
motion is not specified, cosmic rays should travel 
freely in the undefined containment volume. But 
this free motion, a silent element of any variants of 
the Leaky Box Model, implies that a physical, real 
process, at some boundary of the disc, reverses the 
ion motion (reflection). Without the reflection at 
some boundary of the disc the grammage cannot 
accumulate to high values, because the free traver- 
sal of the disc without reflection entails a gram- 
mage of a few milligrams per cm 2 , some 4 orders 
of magnitude below the standard 10 g/cm 2 . What 
is the physical mechanism accomplishing this op- 
eration? To date (2007), it remains unknown. 

Notice further that the matter density in the disc, 
adopted in Leaky Box Models for intrinsic cal- 
culation procedure to determine nuclear spallation 
rates, is also inconsistent with the matter den- 
sity necessary to determine residence time of cos- 
mic rays using radioactive clock measurements 
( 10 Be/Be and others). A tangible sign of this em- 
barassing feature is that the mean gas density in 
the disc turns out to be in the range 0.25 - 0.35 
atom/cm 3 (see, for example, [17] for the data and 
[6] for the calculation procedure), a factor 3-4 be- 
low the average observed value of 1 atom/cm 3 . 

A paradox is encountered by extrapolating the res- 
idence time at energies above 10 17 adopting Gr~ 5 
with (5=0.65 [12]. The high energy galactic sources 
would populate a small volume in the disc, concen- 
tric to the Earth, and they all would reside in the 
solar system at energies above 10 23 eV. 
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